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Invariant mass distributions of the hadronic decay products from resonances formed in relativistic 
heavy ion collision (RHIC) experiments are investigated with a view to disentangle the effects of 
thermal motion and the phase space of decay products from those of intrinsic changes in the struc- 
ture of resonances at the freeze-out conditions. Analytic results of peak mass shifts for the cases of 
both equal and unequal mass decay products are derived. The shift is expressed in terms of the peak 
mass and width of the vacuum or medium-modified spectral functions and temperature. Examples 
of expected shifts in meson (e.g., p, u), and a) and baryon (e.g., A) resonances that are helpful to 
interpret recent RHIC measurements at BNL are provided. Although significant downward mass 
shifts are caused by widened widths of the p— meson in medium, a downward shift of at least 50 
MeV in its intrinsic mass is required to account for the reported downward shift of 60-70 MeV in 
the peak of the p— invariant mass distribution. An observed downward shift from the vacuum peak 
value of the A distinctively signals a significant downward shift in its intrinsic peak mass, since 
unlike for the p— meson, phase space functions produce an upward shift for the A— isobar. 

PACS numbers: 25.75-q, 25.75.Dw, 13.25-k, 13.30.Eg 



I. INTRODUCTION 

Decay products from unstable resonances whose life- 
times are smaller than the duration of ~ 10 — 15 fm/c 
of a relativistic heavy ion collision (RHIC) give valuable 
insights into (a) the ambient conditions of strongly in- 
teracting matter (such as its temperature and chemi- 
cal potential if matter is equilibrated), and (b) the in- 
fluence of the hot and dense medium on their intrinsic 
properties (such as their masses and widths). Through 
their electromagnetic decays, vector meson resonances 
like p, w, di, etc... enable a direct view into hot hadronic 
matter, since dileptons and photons escape without fur- 
ther interactions with matter (for recent accounts, see 
Refs. [1-3]). In contrast, hadronic decays which occur 
when the medium is dense and strongly interacting con- 
tribute to thermalize matter, since the decay products 
undergo substantial rescattering before escaping. Hence 
such hadronic decay products cannot be used to infer the 
properties of the parent resonances. 

An intriguing exception is the case in which resonance 
decays occur close to the thermal freeze-out of hadrons. 
If the decay products do not further interact with the 
hadronic environment, and their total four- momentum 
q can be measured (this requires that the particles in 
the final state can be identified and their momenta mea- 
sured) , a resonance bump should appear in the invariant 
mass plots in M = \/s = \fq* . Such reconstructed reso- 
nances constitute spectroscopic probes of the freeze-out 
stage of RHIC experiments. Where possible, a compar- 
ison of invariant mass plots from e + e~ or pp collisions, 
with those from nucleus-nucleus collisions at the same 
cm. energy enables us to establish the extent to which 
medium effects leave their imprint on resonances. As we 
demonstrate, shifts in the peak masses for the cases of 



p— meson and A— isobar decays are distinctively different 
at the expected freeze-out conditions of RHIC collisions. 

Thanks to the growing multiplicities of final state 
hadrons in RHIC experiments, detailed spectral infor- 
mation is beginning to become available. Examples from 
Au+Au collisions at BNL include the reconstruction of 
(a) i\~*(892)° mesons [4] and 0- mesons [5] at cm. energy 
= 130 GeV, (b) p-decays [6] at ^ii^ = 200 GeV 
(preliminary results), and (c) A, S~,and fl~ decays [7] 
at y / SNN = 130 GeV. Among the results available to date, 
the 60 — 70 MeV downward shift of the p— meson peak 
from its vacuum value of 776 MeV is interesting both be- 
cause such a shift, if confirmed, implies the presence of 
considerable hadronic interactions at freeze-out and be- 
cause it provides the opportunity to pin down in-medium 
properties of mesons that are susceptible to effects of par- 
tial chiral symmetry restoration. 

Our objective in this work is to point out how to ben- 
efit from such detailed information. Toward this goal, 
we study invariant mass distributions from the hadronic 
decay of unstable resonances by exploring the combined 
effects of phase space functions associated with thermal 
motion and the decay products, and possible intrinsic 
changes in the properties of resonances in a hadronic 
heat bath. While phase space modifications can be eval- 
uated without approximation, those induced by intrinsic 
changes of resonance structure require model dependent 
considerations. Ongoing debates about medium modifi- 
cations of spectral functions can be found, for example, 
in Refs. [8-13]. In order to highlight the basic trends of 
in-medium effects on resonance properties, we utilize re- 
sults of model calculations in vogue in conjunction with 
the inevitable phase space functions to show how the in- 
variant mass distributions of resonance decay products 
are modified from their vacuum counterparts. 
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II. INVARIANT MASS DISTRIBUTIONS FROM 
THE DECAY OF UNSTABLE RESONANCES 

In thermalizcd matter, the number of events in which 
a resonance R is produced and subsequently decays into 
a final state / is given by (see, for example, Ref. [14]) 



dN f 



dN R (-Im n)(M fl /7r) 



d i xd i q d 3 xd 3 q (M 2 - M 2 ) 2 + (Im II) 



ouvac 
2 11 A-/ • 



(1) 



Above, q is the total four-momentum of the decay prod- 
ucts, M 2 = s = q 2 is the invariant mass squared, 
dNfi/d 3 xd 3 q is the thermal phase space (Bose-Einstein 
or Fermi-Dirac) distribution of the off-shell unstable res- 
onance particle, Mr is the mass of the resonance at its 
peak, and r™f^ is the decay width in vacuum. The 
proper self-energy II describes the influence of the heat 
bath on the properties of the resonance; explicitly, 



Im n = -M R T 



tot ■ 



M 2 = {M v R ac f + Re n . (2) 



In general, the resonance mass Mr and the width T tot 
depend on the temperature T, chemical potential p, and 
interactions between the various constituents of the heat 
bath. Additionally, both Mr and T to t depend on the mo- 
mentum of the resonance [8]. Eq. (1) applies strictly to 
the case in which the hadronic decay products do not un- 
dergo further interactions within the medium and freely 
escape to the detectors. Its applicability is therefore con- 
tingent on rescattering effects being small at the decou- 
pling or freeze-out conditions. 

The thermal distribution of a resonance particle is 



(2J+ 1) 
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where J denotes the spin of the resonance, = T is 
the temperature, a = 1 for bosons and —1 for fermions, 
and E = Mp cosh y is the energy of the resonance with 
rapidity y and transverse mass Mp = \J M 2 + q\ , in the 
fluid rest frame. 

In what follows, we utilize Eqs. (1), (2), and (3) to 
study the invariant mass distribution of the decay prod- 
ucts stemming from resonances in a mixture of thermal- 
ized hadrons. We begin by first considering p — > -kit 
decays. In vacuum, we take the decay width of the p 
(the decay occurs with an orbital angular momentum of 
I = 1) to be 
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for pion cm. momenta k not too far in excess of ko, 
the value at the peak of the resonance. The width grows 
as fc( 2 ^ +1 ) due to the penetration factor associated with 
a partial wave of orbital angular momentum I. For 
k ko, the monotonically increasing fc— dependence in 



this equation should, however, be tamed as noted e.g., 
in Ref. [15]. Alternative forms for have been used 

in the literature (see, for example, Ref. [16]). Since our 
primary interest is limited to the peak region, Eq. (4) 
provides an adequate representation insofar as these dif- 
ferent parametrizations all agree near the Breit-Wigncr 
peak. The invariant mass distribution is obtained from 
Eq. (1) by integrating out the longitudinal and transverse 
momenta of the resonance, i.e., 
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Equations (1) through (5) enable us to explore the com- 
bined effects of intrinsic changes in the properties of the 
resonance in a heat bath (given by the real and imaginary 
parts of the self-energy II) and phase space functions as- 
sociated with thermal motion and decay products. In 
order to highlight the effects of each, we consider in turn 
the cases in which the resonance mass and width are (a) 
not affected by the heat bath, and (b) influenced by in- 
teractions with the other constituents in the heat bath. 



A. Shifts induced by the phase spaces of thermal 
motion and decay products 

In this case, Re II is set to zero (i.e., no intrinsic mass 
shift) and Tr is taken at its vacuum value. This allows 
us to cast Eq. (5) as 



dNf 



2MM R T R T V ^ 



r 7(M, p, T) ; 



d A xdM (M 2 - M 2 R ) 2 + (MrTr) 2 ' 
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I(M, p, T) = ^3 J dyq T dq T e(E _^ )/T _ ~ ■ (6) 

For M/T 3> 1, it is advantageous to express I above 
(after performing integrations over qx and y) as 
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where the +(— ) sign refers to bosons (fermions) and Ki 
is the modified Bessel function. For invariant masses 
M ~ Mr, where Mr ~ 776 MeV, and temperatures 
in the range 100 < T/MeV < 175, we can retain only the 
j = 1 term in the sum above and employ the asymptotic 
expansion K v (x) ~ (-K/2X) 1 / 2 exp(— x) in order to gain a 
qualitative insight into the shifts caused by phase space 
alone. Collecting the M— dependences from Eqs. (5), (6), 
and (7) together, we arrive at 
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where the factors omitted, including those that depend 
on Mr and Tr, contribute only to the normalization. 
The result in Eq. (8) has the appealing physical interpre- 
tation that the invariant mass distribution is given by the 
simple product "Brcit-Wigncr probability of resonance 
formation" x "Available phase space for decay prod- 
ucts" x "Density of Boltzmann-like 'particles' of mass 
M ." (Such simplicity is lost when effects of interactions 
and Bose-Einstein enhancement are included.) Note that 
while the Boltzmann factor exp(— M/T) tends to pull the 
peak position of the Breit-Wigner at Mr to lower values 
of M (the lower the T, the larger is the downward shift), 
the M 4 5 power law dependence arising from the phase 
space of decay products and the pre-exponcntial factor in 
the density tend to pull the peak position to M > Mr. 
The net result from Eq. (6) is shown in Fig. 1, where we 
have used Mr = 775.9 MeV and T R = 147.9 MeV in- 
ferred from p— production in e+e~ collisions [17] which 
are free of uncertainties associated with p— production in 
hadronic collisions * . The results shown refer to the case 
of /i = 0; nonzero values of p, mainly affect the normal- 
ization of the distributions, but not their peak positions 
or shapes. It is clear that the combined effects of thermal 
and decay product phase spaces are small, and at best 
produce a downward shift of about 10 MeV at the lowest 
temperature in the range 100 < T/MeV < 175. 



*The resonant part of the p— mass spectrum in hadronic col- 
lisions is well fit by a Breit-Wigner form with Mr(p) = 766.5 
MeV and T R (p) = 150.2 MeV [17]. Relative to the spectrum 
observed in e + e~~ collisions, there is ~ 10 MeV downward 
shift in the peak mass, which is commonly attributed to ef- 
fects of final state interactions. While it is tempting to use 
this downward shift to advantage, it is, as yet, unclear if the 
effects of final state interactions in heavy-ion collisions are 
the same as in elementary hadronic collisions. Use of such an 
intrinsic shift is, however, straightforward in the theoretical 
formalism outlined in Sees. II through IV. 
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FIG. 1. p— meson invariant mass distributions in a heat 
bath at the indicated temperatures T. Results exhibit effects 
of thermal and decay product phase spaces alone. The dashed 
line shows the peak p— mass in vacuum inferred from e + e~ 
collisions (see inset for values of Mr and Tr). 

The peak position is determined by locating the value 
of M at which Eq. (6) exhibits a maximum. Setting the 
derivative of Eq. (6) with respect to M to zero, we solve 
for x from 



4x 2 (x 2 - y 2 



+ 



z£°liKo(?V) (2l+l)x 2 



(x 2 -y 2 ) 2 + (yz) 2 £~ iKxO* 



1, 
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where we have introduced the abbreviations x = M/T, 
y = Mr/T, y„ = mn/T and z = Tr/T (for the p—inr de- 
cay, £ = 1). For the temperatures chosen in Fig. 1, Table 
1 gives results for the peak positions M max , shifts from 
the vacuum value AM = M max — Mr, average masses 
(M), and variances ctm = \J (M 2 ) — (M) 2 , where the 
symbol (• • •) denotes averages taken over the invariant 
mass distribution. Together with M max , the quantities 
(M) and <jm provide a rough measure of the shape of the 
distribution. As noted earlier, the shift of the peak mass 
is largest at the smallest temperature, since in this case 
the steep Boltzmann spectrum provides a strong bias to- 
wards low masses. This leads to the intriguing possibility 
that a measured shift of the p-mass can be used to delimit 
the freeze-out or decoupling temperature of the disassem- 
bling hadronic fireball, provided medium-induced effects 
on the p-spectral function and rescattering of the decay 
pions are small or under full theoretical control. 

The opposite pulls of the falling exponential (due to the 
Boltzmann factor) and the rising power law (due to the 
phase space of decay products and the pre-exponential 
factor in the non-degenerate density) nearly cancel at 
T ~ 160 MeV, leading to AM ~ 0. At the highest 
temperature considered (T ~ 175 MeV is thought to be 
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the chemical freeze-out temperature at RHIC [18]), the 
power law rise wins over the exponential fall inasmuch as 
AM > 0. It is worthwhile remembering, however, that 
at these temperatures, effects of interactions cannot be 
ignored, and that the results shown in Fig. 1 and Table 1 
refer to the effects of phase space alone. Furthermore, as 
T — > T c , the basic structure of a resonance itself would be 
considerably altered from its vacuum form due to effects 
of deconfinement and/or chiral symmetry restoration. 



Note that the second and third terms in Eq. (11) give sub- 
stantial contributions, which attests to the significance 
of the power law behaviors of the available phase space 
(see Eq. (4)) and the pre-exponcntial term in the den- 
sity of particles of mass M. For example, the first term 
A M = -r 2 /(8T) in Eq. (12) yields a shift of 

AM = -27.3 (-15.6) MeV for T = 100 (175) MeV , 

(13) 



T 
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150 
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717.7 
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Table 1: Basic characteristics of the p— invariant mass dis- 
tributions in Fig. 1. T is the temperature, M max is the peak 
position, AM = M m ax — Mr is the shift from the vacuum 
value, (M) is the average mass, and om = \J (M 2 ) — (M) 2 is 
the variance of the distribution (all values are in MeV). The 
vacuum Breit-Wigner peak mass and width are Mr = M p — 
775.9 MeV and T R = V P = 147.9 MeV, respectively. 



which is considerably larger in magnitude than -9.8(1.6) 
MeV given by considering the other terms. The corre- 
sponding shift from Eq. (12) is -9.2(1.9) MeV. These 
shifts are to be contrasted with the full numerical so- 
lutions from Eq. (9): 

AM = -10.6 (1.5 MeV) for T = 100 (175) MeV . 

(14) 

Observe that Eq. (11) implies that AM = for 
Mr 



160 MeV . 



(15) 



B. Analytic results to leading order 

To better understand the dependence of the shift on 
the decoupling temperature, peak mass, and width of the 
resonance, we solve Eq. (9) to linear order in the shift 
Sx = x — y, which is at best of order 10% (see Table 1). 
Keeping only the first term in the sums in Eq. (9) (the 
Boltzmann approximation is adequate around the peak 
mass), we find 
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where the Bessel functions are evaluated at y — Mr/T. 
For y^ 1 = T/Mr <C 1, we can simplify the above expres- 
sion to obtain 



AM 



T 
Mr 



8T 2 /r 2 + T/M R 



(11) 



In obtaining this result, use was made of the asymptotic 
expansion K (y)/Ki(y) ~ 1 — l/(2y) + ••• . Further 
linearization of the fraction in Eq. (11) in terms of T /Mr 
results in 
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8T 
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in excellent agreement with that obtained from Eq. (9). 
A comparison of M max versus temperature from the dif- 
ferent approximations is shown in Fig. 2. We have veri- 
fied that effects of Bose statistics are not very significant 
at these temperatures. 
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FIG. 2. Peak positions of the invariant mass distribution 
functions in Fig. 1 versus temperature. The solid curve refers 
to the exact numerical results from Eq. (9), while the dashed 
and dash-dotted curves are for the approximations in Eqs. 
(11) and (12), respectively. The dotted curve refers to the 
zeroth order (in T/Mr) result A Mo = -r 2 /8T. 
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C. Shifts induced by in- medium interactions plus 
thermal and decay product phase spaces 

From Eqs. (11) and (12), it is evident that the width 
of a resonance strongly influences the magnitude of the 
peak mass shift. Thus, a broad resonance such as the 
/o(600) will exhibit a larger shift than sharp resonances 
such as the K* or <f>. In-mcdium interactions with the 
other constituents of the heat bath, namely mesons such 
as 7r, K, .., u>, <fi, ai, ... and baryons such as A, TV*..., will 
likely modify the p— spectral function from its vacuum 
form. The observed distributions are, however, convolu- 
tions of such intrinsic changes (pressure broadening and 
interaction-induced mass shift) with the thermal phase 
space and the phase space for the corresponding decay 
products. These intrinsic changes are given by the self- 
energy II (see Eq. (2)), with Re II giving the intrinsic 
mass shift and Im II giving the in-medium width. 

Before turning to detailed model calculations of II, it is 
useful to explore the combined effects of schematic intrin- 
sic changes (such as a broadened, but unshifted resonance 
or a shifted, but not broadened resonance; the realistic 
case would include a combination of both) and the ther- 
mal and decay product phase space functions. In Fig. 3, 
we show the invariant mass distribution of the p-meson 
with a fixed peak mass of 775.9 MeV at a temperature 
of 100 MeV for different widths (T R = 100 to 200 MeV 
in steps of 25 MeV) . It is evident that an intrinsic broad- 
ening of the resonance induces a considerable shift in the 
peak mass; the larger the broadening, the larger is the 
shift. 
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FIG. 3. Invariant mass distributions of the p meson in a 
hadronic heat bath (at a temperature of 100 MeV) for vary- 
ing intrinsic widths. The dashed line indicates the unshifted 
rho-mass Mr = M p = 775.9 MeV. 

It is noteworthy that for low temperatures and large 
intrinsic widths, the shape of the invariant distribution 
is considerably distorted from its symmetric form in vac- 



uum. Consequently, the definition of the "width" as nor- 
mally used for symmetric shapes loses its utility as shown 
by the standard deviations um in Table 2. The locations 
of the maxima and the mass-shifts are also listed in Table 
2 for the various curves shown in Fig. 3. 
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Table 2: Shifts (from Eq. (9)) of the peak mass induced by 
varying intrinsic widths, Tr, of the resonance. The Breit- 
Wigner peak mass is at Mr — M p — 775.9 MeV. Also shown 
are the average mass and the standard deviation of the dis- 
tribution. All values are in MeV. 



We turn now to the case in which intrinsic changes in 
the peak mass alone are considered. For simplicity, and 
because the dependence on the momentum of the reso- 
nance is small for low temperatures and baryon number 
densities at freeze-out (see Fig. 8 of Ref. [8]), we use con- 
stant peak mass shifts for illustrative purposes. The re- 
sults in Fig. 4 and Table 3 confirm our expectations from 
Eqs. (11) and (12) that the mass shift is largely indepen- 
dent of the peak mass. For the case considered, T = 100 
MeV, the mass shift relative to the peak mass varies from 
1.4% for M R = 836 MeV to 1.1% for M R = 716 MeV. 
Thus, relative to the assumed peak mass, shifts induced 
by intrinsic mass shifts are smaller than those caused by 
changes in the width. 




100 
, R (MeV) 

FIG. 4. Invariant mass distributions of the p— meson with 
a width r_R = r p = 147.9 MeV and assumed intrinsic peak 
masses that vary from 716 to 836 MeV. For each curve, Mr 
is subtracted from the horizontal scale. 
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Table 3: Peak mass M max and shift AM = M max — M R (from 
Eq. (9)) versus the Breit-Wigner maximum Mr for a fixed 
width T R = T p = 147.9 MeV (all numbers are in MeV). 
Also shown is the the standard deviation om of the invariant 
distribution. 

The results shown in Fig. 4 and Table 3 refer to peak 
mass shifts caused by the phase space functions on top 
of the assumed intrinsic mass shifts in the range —60 < 
AM < 60 MeV t. Compared to the vacuum peak mass, 
the total shift is given by 



AM tot = {M R - 775.9) + AM . 
Thus, for example, 



(16) 



AM tot = -68.1 MeV for M R = 716 MeV , 

= +48.3 MeV for M R = 836 MeV . (17) 

It is evident that sizeable intrinsic mass shifts are re- 
quired to generate sizeable peak mass shifts in the in- 
variant mass distributions. 



III. CONTRIBUTIONS TO THE nn FINAL STATE 
FROM THE ^(782) AND F (600) MESONS 

In the peak region of the p— invariant mass spectrum, 
contributions to the itit final state include those that arise 
from the f (600) (hereafter a) and u> mesons that are also 
likely modified by medium effects. 



A. lu — > ir + ir decays 

The decay of the lu is distinctive, chiefly because of its 
narrow total width T = 8.44 MeV in vacuum. Although 



the dominant decay channel is the lu 



decay 



(iyr = 89.1 ± 0.7%), a small fraction (iyr = 1.70 ± 
0.28%) of the decays occur via the lu — > n + ir~ channel. 
For the latter case, Eqs. (9) through (12) show that the 
influence of the phase space functions on the peak mass 
shift of the u> is negligibly small. Unless large intrinsic 



shifts are caused by the medium, the shift in the lu— peak 
mass in the invariant distribution of its tt + it~ final state 
is likely to be very small. Calculations to date [8,9] have 
found only modest changes in Re n w and Im II W at the 
freeze-out temperatures of (100-125) MeV. 



B. a — > 7r + 7r decays 

Pinning down the contributions from the a decays 
is more difficult, especially since its vacuum width is 
rather large (600 — 1000 MeV) and its peak position 
(400— 1200 MeV) in vacuum is somewhat uncertain [17]. 
To encompass the possible variations, we first present 
results of two illustrative calculations; one in which we 
take M R {a) = 800 MeV and T R (a) = 800 MeV (this 
case refers to the average peak mass and the average 
width, respectively), and the other in which we employ 
M R (a) = 600 MeV and T r (ct) = 400 MeV (these val- 
ues represent averages from the D— and r— decays [17]). 
The vacuum decay width for the a — > 7T7t is 



r^ /( .) = r fl (.)(A 



2^+1 



M 2 -4m 2 \ 1/2 
M^a)-4ml 



(18) 



as, in this case, the orbital angular momentum I = 0. 
Utilizing this in Eq. (6) with spin J = 0, results of in- 
variant mass distributions from the a — > -kit decay are 
shown by the dash-dotted lines in Fig. 5 for tempera- 
tures expected near freeze-out. As is clear from the solid 
curves in this figure, which show the summed contribu- 
tions from the p — > 7T7T (dashed curves) and a — > itit de- 
cays, the maximum near the peak-region of the p— mass 
spectrum is barely influenced by the small and exponen- 
tially falling contributions from the a — > itit decays. The 
contributions from the a — > itit decays dominate, how- 
ever, for masses well below the p— peak. 

By setting I = in Eq. (9), one can straightforwardly 
obtain a very rough estimate of the peak mass shift in 
the invariant mass distribution from the a — > tttt decay. 
The result to linear order in T/M R is 



AM 



T 

Mr 



8 T 2 /r 2 + T/M R 



(19) 



This result, however, provides a poor approximation to 
the numerical results from Eq. (9) chiefly because the 
shifts are large. 



f Upward intrinsic mass shifts are also considered to encom- 
pass the expectations of Ref. [12]. 
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FIG. 5. Invariant mass distributions from the decays of 
the p (dashed curves) and a (dash-dotted curves) mesons. 
The solid curves show the sum of the contributions. The 
Breit-Wigner maximum masses and widths are shown in the 
insets. 

Table 4 lists the peak masses M max and shifts AM 
calculated numerically from Eq. (9) with i = for the 
cases M R (T R ) 600 (400) and 800 (800) MeV, respectively. 
Notice that the shifts in the invariant mass distributions 
from the a— decay in medium depend strongly on the pa- 
rameters used to describe the shape of the a— resonance 
in vacuum. 





T 


-^max 


AM 


Mr 


= 600 


125 


527.2 


-72.8 


Tr 


= 400 


100 


458.5 


-141.5 


Mr 


= 800 


125 


450.5 


-349.5 


Tr 


= 800 


100 


391.0 


-409.0 



Table 4: Peak mass M max and shift AM = M max - Mr (from 
Eq. (9)) versus temperature T for the indicated values of the 
Breit-Wigner maximum Mr and width Fr (all numbers are 
in MeV). 



C. In-medium shifts of the a— meson 

The intrinsic mass shift of the fragile cr— meson is inti- 
mately related with chiral symmetry restoration (cf. [19] 
and references therein) . With increasing temperature (or 
baryon density), Mr{<j) is expected to steadily decrease 
and approach the pion mass m w as chiral symmetry is 
fully restored at T c ~ 170 ± 10 MeV (the precise number 
is still a matter of intense debate!). Being a Goldstone 
boson, niir remains nearly constant at its vacuum value 
even as T c is approached. Thus, at the freeze-out temper- 
atures T ~ 100 — 125 MeV and rather small baryon den- 
sities under consideration here, modest downward shifts 
of the cr— peak mass from its vacuum value are to be 
expected. For definiteness, we choose an optimistic in- 
trinsic mass shift of 100 MeV relative to its vacuum value 
of 600 MeV based on the model calculations reported in 
[19]. At first, we take T R {a) = 400 MeV. In this case, the 
shifts relative to Mr = 500 MeV are AM = -73(-37) 
MeV for T = 100(125) MeV, respectively. These shifts, 
relative to the vacuum peak mass of 600 MeV, amount 
to the total shifts 

AM tot = -173 MeV for T = 100 MeV , 

= -137 MeV for T = 125 MeV . (20) 

In addition to its mass, the cr— decay width Tr{(j) also 
decreases with increasing temperature. For illustration, 
we employ Tr(<t) = 300 MeV. This results in shifts 
AM = -46.9(-23.2) MeV for T = 100(125) MeV. The 
corresponding total shifts are 

AM tot = -146.9 MeV for T = 100 MeV , 

= -123.2 MeV for T = 125 MeV . (21) 

In all cases considered, the a — > 7T7t contributions at 
the peak p— mass are given predominantly by the expo- 
nentially falling Boltzmann factors and provide small cor- 
rections to those from p — ► 7T7t decays. The left-wings of 
the p— distributions are, however, destroyed insofar as the 
summed invariant distributions are relatively flat over a 
wide range of invariant masses M < 4.5™^. Note, how- 
ever, that the if q (498) — ► 7T7t decays that essentially oc- 
cur in vacuum exhibit a prominent peak at ~ 500 MeV 
and that the a — > 7T7t decays constitute a background in 
the vicinity of this peak. 

IV. SHIFTS OF BARYON RESONANCES 

Most of the secondary particles produced in RHIC col- 
lisions are mesons. However, many baryons and baryon 
resonances including their anti-particles are also pro- 
duced. Invariant distributions of hadronic decay prod- 
ucts from baryon resonances are best exemplified by 
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considering the decay of the A(1232)-resonance. In- 
medium properties of the A have been studied earlier 
through photo-absorption and hadronic reactions on nu- 
clei [20-22] , but has received relatively little attention at 
temperatures characteristic of RHIC collisions. 



A. Effects of thermal motion and the phase space of 
unequal mass decay products 

As in the case of p — > irir decay, the A — > nir decay, 
where n is a nucleon, occurs with orbital angular momen- 
tum I = 1. Thus, the invariant mass distribution can be 
calculated following the procedure outlined in Sec. II, 
but with 



pvac 



k J 



2i+l 



(22) 



where k is the cm. momentum of either the pion or 
nucleon, i.e., 



_ 1 / (M 2 - Ml){M 2 - Af2)\ 1/2 
k ~ 2 \ M 2 



(23) 



where M± = M n ± and ko is the cm. momentum 
with M = Mr. (Near the peak region, the functional 
form in Eq. (22) matches that advocated in Ref. [22] in 
which a more complete discussion of Y^^j(A) can be 
found.) In this case, the maximum in the invariant mass 
distribution is obtained by solving for x from 



T 


AM(Eq. (24)) 


AM(Eq. (25)) 


175 


12.6 


13.8 


150 


11.0 


12.0 


125 


8.7 


9.6 


100 


5.4 


6.1 



Table 5: Comparison of peak mass shifts AM = M max — Mr 
from Eq. (24)) and Eq. (25)) versus temperature T for the 
A" — ► pir~ decay. The Breit-Wigner maximum Mr = 1232 
MeV and Y R = 120 MeV. The shifts are in MeV. 



The temperature at which the shift AM = is given 



by 



Mr 

|+3Q 



74.5 MeV. 



(27) 



This temperature, being considerably lower than the 
freeze-out temperature 100 < T//MeV < 125, implies 
that an observed downward shift from the vacuum peak 
value of 1232 MeV from A— decay distinctively signals 
an intrinsic downward shift due to in-medium effects 
By the same token, an observed upward shift must be 
attributed to the predominance of the phase space func- 
tion associated with the decay products, although the 
left-wings of higher mass N* resonances can also con- 
tribute in determining the peak position. We turn now 
to effects of in-medium mass and width modifications on 
the peak mass shift. 



Ax 2 (x 2 - y 2 ) 
(x 2 — y 2 ) 2 + (y z) 2 

3(x 4 



x T,f=i K oU x ) 

£°°=i7 K i(i*) 
- vlv-) _ , 



{x 2 - y%){x 2 - yt) 



(24) 



where x = M/T, y = M R /T, y± = M±/T, and z = 
Yr/T. Since for masses near the peak region, M/T 3> 1, 
Eq. (24) can be solved for the linear shift Sx = x — y as 
in Sec. II. B by considering only the j = 1 terms in the 
second term (the Boltzmann limit). To linear order in 
T/Mr, the shift in the peak mass is given by 



r 2 

AM ~ 

8T 



3 Y 2 \ T 

^ { 2 +3Q+ xF 2 )m-r 



where 



Q = 



M R - (M+M_) 2 
{M 2 R - M 2 ){M 2 R - M 2 ) 



(25) 



(26) 



Note that Eq. (25) reduces to Eq. (12) in the limit of 
identical final state particles. The extent to which the 
linear shift in Eq. (25) reproduces the numerical results 
from Eq. (24) is shown in Table 5. 



B. In-medium shifts of the A isobar 

In photonuclear reactions, cross-sections per nucleon 
for nuclei ranging from 9 Be to 238 U when compared to 
the average single nucleon total photoabsorption cross- 
section clearly show that the A— resonance appears at 
an energy close to that of the free A (1232), but spread 
over a wider energy range [23]. In tt— nucleus reactions, 
however, in addition to a widened width, the resonance 
energy is also shifted downward by about 20-30 MeV. An 
illuminating discussion of the various effects contributing 
to in-medium shifts at zero temperature can be found 
in [22]. The total width arises from a combination of 
competing effects: 



r(A) 



pfree + pel 



^abs _|_ jpPauli 



(28) 



*The precise value of T z depends on the peak mass and width 
of the A— resonance. Intrinsic shifts of these quantities con- 
sidered in Sec. IV.B result in values of T z ~ 75 ± 25 MeV. 
Thus, the main conclusion that T z < Tf remains valid in the 
presence of medium modifications of the width and peak mass 
of the A— resonance. 
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The dominant contribution arises from the so-called elas- 
tic broadening of the width, r cl , given by the imaginary 
part of the matrix element ((Ah)p\Vx\(Ah) a ), where 
h stands for a nucleon hole and denotes the pion- 
mediatcd interaction between baryons 1 and 2. Cou- 
pling to absorptive channels (typically one-quarter to 
one-third of the total cross-section in the resonance re- 
gion comes from this source) corresponds to the coupling 
of A— hole states to multiple nucleon-hole (nNnh) states 
and increases the width. In contrast, the width is re- 
duced by the Pauli exclusion principle, since in the decay 
(Ah) — ► n(Nh), occupied nucleon states are inaccessible. 
The net effect, however, is a net increase in the width 
of (Ah) states. For low baryon densities and for tem- 
peratures of interest here, however, the dominant contri- 
butions arise from T cl and r abs , since Pauli effects are 
small for non-degenerate baryons. Additional modifica- 
tions from coupling to higher-lying N*- resonances remain 
largely unexplored. 

The width modifications are accompanied by disper- 
sive energy shifts (for nuclei, these amount to a few tens 
of MeV). Such shifts are due to the combined effects 
of A— binding in matter, dispersive shifts related to ab- 
sorption, short-range A— hole correlations, coupling to 
higher-lying iV*-resonances, etc. To date, the individual 
contributions from these various sources have not been 
resolved in experiments. 





T AM AM tot 


Mr = 1192 


125 14.2 -25.9 
100 10.9 -29.1 


Mr = 1212 


125 11.1 -8.9 
100 7.9 -12.2 


Mr = 1252 


125 6.8 26.8 
100 3.5 23.5 


Mr = 1272 


125 5.1 45.1 
100 1.8 41.8 



Table 6: Peak mass shifts AM from Eq. (24) caused by 
in-medium Breit-Wigner peak values of Mr and the vac- 
uum value of r fl = T A = 120 MeV. The values of T 
bracket the likefy temperature at freeze-out. AMtot = (Mr — 
1232 MeV) + AM refers to the totaf shift from the vacuum 
peak vafue. Aff numbers are in MeV. 





T AM 


Tr = 100 


125 6.3 
100 3.9 


Tr = 140 


125 11.5 
100 7.1 



Tabfe 7: Peak mass shifts AM from Eq. (24) due to in- 
medium Breit-Wigner widths Tr and the vacuum peak vafue 



of Mr = M A = 1232 MeV. The vafues of T bracket the likefy 
temperature at freeze-out. Alf numbers are in MeV. 

Tables 6 and 7 show results of shifts AM and AMtot = 
(Mr - 1232 MeV) + AM calculated from Eq. (24) with 
representative in-medium shifts. These results show how 
medium-induced modifications are reflected in the invari- 
ant distribution of pir final states from A— decays alone. 
Viewed together with the results in Table 5, the results 
in Tables 6 and 7 imply that a downward shift in the 
peak value of the A— resonance requires a substantial 
dispersive energy shift in the medium. Building upon 
the known wNA physics, an in-depth finite temperature 
analysis that includes coupling to higher N* resonances 
is necessary. 

V. DISCUSSION 

A. In-medium spectral functions 

In the context of lepton-pair (e.g., e+e~) and photon 
production from the dense stages of heavy-ion collisions, 
considerable amount of work has been done on the in- 
medium spectral functions of vector mesons (see, e.g., 
Refs. [2,3] for exhaustive references). The viewpoints 
that have emerged can be roughly classified into two 
broad categories: 

(1) The masses of the vector mesons, particularly those 
of the p and u), decrease substantially with increasing 
temperature T and net baryon density rig. This view, 
espoused in Ref. [13] and shared to varying degrees in 
other works [19,24], is based essentially on the decreasing 
behavior of the (qq) condensate which in turn is related to 
the vector meson masses as a function of T and ub- For 
additional dependences on strong interaction couplings, 
see Refs. [13,19,24]. 

(2) The peaks of the spectral densities are little shifted 
from their vacuum positions, but the widths are consid- 
erably increased due to medium effects that include col- 
lisional broadening [2,3,8]. This view stems from many- 
body calculations that consider the many hadron states 
(for an extensive list, see Sec. II of Ref. [8], and par- 
ticularly, Table 1 therein) to which the vector mesons 
couple in a heat bath. By relating the self-energy to the 
real and imaginary parts of the forward scattering ampli- 
tude, Ref. [8] exploits experimental data where possible 
to infer the peak positions and widths in medium. 

In both scenarios, the largest medium-induced effects, 
whether it is a mass shift or a widened width, occur for 
T -> T c ~ 170 ± 10 MeV and for n B > n , where n ~ 
0.16 fm -3 is the nuclear equilibrium density. However, 
near freeze-out T ~ (100-125) MeV < T c and n B < n . 
Since our focus here is on the properties of resonances as 
reflected by their hadronic decay products at the dilute 
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stages of freeze-out, medium-induced effects are expected 
to be considerably smaller than those for lepton-pair and 
photon production which occurs dominantly at the dense 
stages of the collision. 

In computing shifts of peaks in the invariant mass dis- 
tributions in this work, we have utilized the basic trends 
found in the intrinsic shifts of both the peak masses and 
widths in both of the above approaches. Specifically, the 
effects of widened widths and both downward and up- 
ward intrinsic mass shifts on the invariant mass distribu- 
tions were investigated. For temperatures at freeze-out 
and negligibly small baryon densities under considera- 
tion, the momentum dependences of the intrinsic widths 
are rather weak [8] being significant only for momenta 
q > Qav — VSMT. These weak momentum dependences 
are washed out upon convolution with thermal and de- 
cay product phase space functions. Thus, our discussion 
has been carried out with momentum-independent in- 
medium widths. For the decay widths, however, the ap- 
propriate orbital angular momentum dependences were 
considered. 



B. Flow— related matters 

Hadrons of various species (n,K,N, etc.,) are known 
to exhibit strong collective expansions (large longitudinal 
and transverse flows have become a part of RHIC life; see 
[25] and references therein) as a result of pressure gradi- 
ents established at a relatively early stage of the collision 
[26] . Prime examples of collective flow include the "blue- 
shift" or flattening of the observed transverse momentum 
spectra [27] , and the azimuthal angular anisotropy or el- 
liptic flow [28], of the various particles. Invariant mass 
spectra, however, do not depend on Lorentz boosts and 
are thus unaffected by collective flow. Contributions from 
the various phase space cells can be determined in their 
rest frames and subsequently summed to give the total 
distribution. This underscores the importance of spec- 
tral function studies that are potentially the only means 
to delineate the intrinsic properties of resonances at both 
the dense and freeze-out stages of high-energy heavy-ion 
collisions. 



C. Relation to other works 

Based on the nonrelativstic expression for the Breit- 
Wigner peak 



AM ~ - 



1 



p(M) 



Mr 



(M-M K ) 2 + r 2 /4 



(29) 



Ref. [29] found a downward energy shift of the maximum 
of the energy spectrum of the p— mesons to be 



3 T 
8T V~ 2lfe, 
17.5 MeV for T 



120 MeV . 



(30) 



This shift differs substantially from the —5.4 MeV shift at 
the same temperature obtained in this work. It is easy to 
verify that shifts to linear order obtained from the use of 
nonrelativistic and relativistic Brcit-Wigner forms both 
lead to the same result in Eq. (10). We therefore surmise 
that the substantial positive shift from the third term in 
the numerator of Eq. (11), 



(AM) 



phase space 



3 1 



Ami 



T 

8T ' V" M\) ' Mfl 
12.2 MeV for T= 120 MeV. (31) 



which arises from the phase space of the decay products, 
was possibly not fully considered in estimating the shift 
in Ref. [29]. 

Our work here partly overlaps that in Ref. [30]. Com- 
plementary information, including mass shifts through 
the Brown-Rho scenario, the role of t— channel reso- 
nances in inducing mass shifts of the p— meson, some 
aspects of kinetic freeze-out, effects of the heat bath etc., 
can be found in Ref. [30]. However, the influence of 
the angular-momentum-dependent decay product phase 
space, crucial for comparison with experimental invari- 
ant mass distributions, was not considered in [30]. Our 
work, in addition to bracketing the effects of different sce- 
narios concerning medium-induced effects, highlights the 
importance of the competition between intrinsic shifts 
and phase space functions related with thermal effects 
and decay products. The large differences in the peak 
mass shifts in invariant mass distributions between our 
work and those of Rcfs. [29,30] stem primarily from the 
fact that the decay product phase space functions were 
not considered in the latter works. 

The results in Sees. II through IV clearly show 
that qualitative differences exist in the influence of the 
M— dependence on the £ = 1 decays of p —> irir and 
A — > ]37r, insofar as the decay product phase space func- 
tions, in conjunction with thermal phase space and ex- 
pected intrinsic shifts, produce qualitatively different re- 
sults. In p — > 7T7T decay, the net result is a downward shift 
from the vacuum value of the p— peak mass at freeze-out 
temperatures. In A — > pir decay, however, larger than 
vacuum peak masses are favored unless large modifica- 
tions in the intrinsic mass and width tip the balance. 
Our work shares the view with that of Ref. [30] that 
substantial in-medium collisional and dispersive energy 
shifts are necessary to produce the substantial mass shift 
in the reported p— meson invariant mass spectrum [6]. 
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VI. SUMMARY AND PRINCIPAL FINDINGS 

A study of the peak positions of resonances and their 
shapes is a time-honored spectroscopic tool that is be- 
ginning to be used to probe the freeze-out conditions at 
RHIC collisions and to establish the properties of reso- 
nances in a hot and dense medium. The recently reported 
60-70 MeV downward shift in the peak of the invari- 
ant mass distribution of pions from the decay of p(776) 
mesons is particularly interesting in this regard, since the 
decay products originate from the dilute freeze-out stage 
of the collision. In this work, we have investigated in 
some detail the interplay of thermal motion, the phase 
space of decay products, and intrinsic changes in the 
structure of the p(776), w(782), /o(600), and A(1232) 
resonances for temperatures and hadron densities ex- 
pected at freeze-out. The effects of thermal and or- 
bital angular momentum-dependent phase space func- 
tions were evaluated without approximation. Analytic 
results of shifts to leading order in T/Mr (T is the tem- 
perature and Mr is mass at the peak of the resonance) 
for the cases of both equal and unequal mass decay prod- 
ucts were derived. Medium-induced intrinsic shifts in the 
width and peak positions were taken from state-of-the- 
art calculations in order to assess the extent to which 
such shifts are required to account for the reported shift 
of the p— meson peak. The contribution to the 2tt final 
state from the /o(600) or sigma meson was assessed. The 
expected shift in the peak of the invariant mass distribu- 
tion of decay of the A(1232) resonance was calculated. 

Among the important points that emerge from calcu- 
lations in which resonance peak masses and widths were 
taken to be at their vacuum values are: 

• The combined effects of thermal and decay product 
phase spaces on the p— peak mass are small, and 
at best produce a downward shift of about 10 MeV 
(from the vacuum value) at the lowest temperature 
in the range 100 < T/McV < 160. Upward shifts 
occur for temperatures T > 160 MeV. 

• The (7 — invariant distribution is considerably dis- 
torted, chiefly because of its large vacuum width. 
Depending strongly on the parameters used to de- 
scribe its mass spectrum in vacuum, shifts of hun- 
dreds of MeV occur in its peak mass in a much 
broadened mass distribution. 

• The shift in the peak mass of the co— resonance is 
negligibly small because of its small vacuum width. 

• Largely due to the growing phase space of the de- 
cay products, the peak of the A— resonance shifts 
upward from its vacuum value in the range 100 < 
Xy/MeV < 120, where Tf is the freeze-out tem- 
perature. Downward shifts are produced only for 
(unlikely) freeze-out temperatures Tf < 75 MeV. 



For freeze-out temperatures in the range 100 < 
Tf/MeV < 125, calculations in which intrinsic changes in 
the resonance widths and masses were incorporated led 
to the following conclusions: 

• Although significant downward mass shifts are 
caused by widened widths of the p— meson, a size- 
able shift in its intrinsic mass, at least of order 50 
MeV, is required to generate a shift of order 60-70 
MeV in the peak of the invariant mass distribution 
of the p— meson. 

• Downward shifts of hundreds of MeV occur in 
the o — invariant distribution with decreasing in- 
medium a masses suggested by chiral symmetry 
restoration. As a result, contributions from the 
a — > tttt decay at the peak p— mass are primarily 
determined by the exponentially decreasing ther- 
mal phase space and provide small corrections to 
those from p — > tttt decays. The summed in- 
variant distribution is relatively flat in the range 
2m- K < M < 4.5m v . In this region, however, the 
A^q(498) — ► tttt decays exhibit a prominent peak at 
<~ 500 MeV and the a — ► tttt decays constitute a 
background to this peak. 

• Unless large intrinsic shifts arc caused by the 
medium, the shift of the lo— peak in the invariant 
distribution of its tt + tt~ final state is very small. 

• An observed downward shift from the vacuum peak 
at 1232 MeV from A— decay distinctively signals 
a downward shift in its intrinsic peak mass due 
to in-medium effects. (Recall that, unlike for the 
p— meson, thermal and decay product phase space 
functions produce an upward shift for the A.) For 
example, at Tf = 100 MeV, the peak mass in the in- 
variant distribution moves down by <~ 10(30) MeV 
for an intrinsic downward mass shift of 20(40) MeV. 

VII. OUTLOOK 

Some points that are naturally raised as a consequence 
of this work and that merit further study include: 

• The reported shift in the p— meson peak is sur- 
prisingly large. Why and how then are the parti- 
cles decoupling insofar as freeze-out represents the 
stage when interactions are at their weakest? The 
answers to these questions hinge on the extent to 
which rescattering effects affect the decay distribu- 
tions at freeze-out (see, e.g., Ref. [29]). 

• A shift as large as 60-70 MeV in the p— mass at 
the dilute stages points to considerably larger shifts 
at the dense stages of the collision. This has 
far-reaching consequences for dilepton and photon 
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emissions, which occur predominantly at hadron 
densities and temperatures that are significantly 
larger than those at freeze-out. Experimental es- 
tablishment of spectroscopic shifts, both at the di- 
lute stage (through hadronic probes) and dense 
stage (through electromagnetic probes) at RHIC 
collisions offers the opportunity to choose between 
scenarios that have distinctive predictions about 
the effects of chiral symmetry restoration on res- 
onances in a hot and dense medium. 

• Spectroscopy of resonance decays, particularly 
those of the p— meson and A-isobar, also offers 
the means to delimit the temperature and baryon 
chemical potential of hadronic matter at the decou- 
pling stage (for more details, see Sees. II and IV). 
This complements our knowledge about freeze-out 
conditions inferred from stable particle ratios. 

• We cannot refrain from mentioning that it would 
be marvelous, albeit difficult, to conduct resonance 
tomography through Hanbury-Brown-Twiss anal- 
yses by utilizing the decay products of resonances 
at RHIC and planned LHC experiments. 
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